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We discuss the finite-size efi'ects on the chiral phase transition in QCD. We employ the Nambu- 
Jona-Lasinio model and calculate the thermodynamic potential in the mean-field approximation. 
Finite-size effects on the thermodynamic potential are taken into account by employing the multiple 
reflection expansion. The critical temperature is lowered and the order of the phase transition is 
changed form first to second as the size of the system of interest is reduced. 
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_ I. INTRODUCTION 

^ . 

["t I , The phase structure of QCD at high temperature T and/or quark chemical potential /i has attracted a great deal 
of interest in cosmology, compact stars and heavy-ion collisions. During the last decades, significant advances have 
been made in our understandings of the phase structure of hot and/or dense quark matter. At the present time, 

■ it is widely accepted that the QCD vacuum undergoes a phase transition (into a chirally symmetric phase or a 
^ \ color-superconducting phase) at sufficiently high T and /i. 

1 The studies of the phase structure of QCD so far have been exclusively devoted to homogeneous quark matter in 
5l5 bulk. However, the QGP phase created in relativistic heavy-ion collisions has a finite size and it is not obvious that 
the size is large enough to allow to apply the thermodynamic limit. If the size of the system of interest is not large 
enough, we need to take account of the deviations from thermodynamic calculations. For instance, the fluctuations 
of order parameter, induced by the finite-size effects, would lower the critical temperature and change the order of 

■ phase transition. 

I There are many works which concern with finite-size effects on hadron physics; lattice QCD calculations 
. 13, I M B. Ip , the finite-size scaling analyses 1^, 0, ^| , plasmon [13] , baryon number fluctuations , and so on 
I ; ]l Jlia lld Tll .18. 19- 20- 21, 22, 23, 24= 2a,!26[. In particular, the finite-size effects on the QCD phase diagram are 
Qh, studied in [23| and j2J|. In the former, however, the discussion is limited to (2-f l)-dimensional systems and in the 
^ ' latter, the phase transition discussed is not related directly with the manifestation of chiral symmetry. 
'T^ I In this paper, we discuss the finite-size effects to the chiral phase transition. For the purpose of illustrating our 
K*" ■ approach, we adopt the two-flavor Nambu-Jona-Lasinio (NJL) model as an effective model of QCD. The finite-size 
k>( ' effects are taken into account by using the multiple reflection expansion (MRE) 113, IH,!!^. Then, we study the size 
j_j dependence of the thermodynamic potential in the mean-field approximation. 

d ' In the MRE approximation, the finite-size effects are included in terms of the modified density of states. The MRE 
has been used to calculate the thermodynamic quantities (such as the energy per baryon and the free energy) of finite 
lumps of quark matter js^l • As far as overall structures are concerned, the results are in good agreement with those 
of the mode-filling calculations with the MIT bag wave functions. We emphasize that our model calculation differs 
from the finite-size scaling analyses (in which finite-size effects act as an external field; more specifically, mass) and 
the studies in a finite box in that the present model is essentially based on the MIT bag wave functions. However, it 
should be noted that the relevance of the MRE approximation is still inconclusive when it is applied to nonperturbative 
calculations [^H^, 33, 34]. 

This paper is organized as follows. In section 2, we apply the MRE approximation to the NJL model and calculate 
the T-/X phase diagram. However, the result is unacceptable because the order of chiral phase transition is changed 
from second to first. In section 3, we reconsider the application of the MRE to the chiral phase transition. Summary 
and concluding remarks are given in section 4. 
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II. SIMPLE APPLICATION OF MRE 



We consider up and down quarks to be massless and start with the foUowing SU(2)i x SU(2)ii; symmetric NJL 
Lagrangian, 



£ = qi^dt^q + G [{qqf + [qil^Tqf 



(1) 



where q denotes the quark field with two flavors [Nj — 2) and three colors {N^ — 3) and G is the coupling constant. 
The Pauli matrices r act in the flavor space. The ultraviolet cutoff Auv = 0.65 GeV and the coupling constant 
G = 5.01 GeV~^ are determined so as to reproduce the pion decay constant /jr — 93 MeV and the chiral condensate 
{qq) = (-250 MeVf in the chiral limit Ullll- 
In the mean-field approximation, the thermodynamic potential per unit volume lo is given by 



m 
4G 



^uv 



dkp{k)i^Ek+'^\n 



1 



1 + g-ZSlfifc-p) 



(2) 



where v = 2NcNf, (3 = 1/T and Ek = ^JW 
The density of states p{}i) is given by 



rn? with m ~ ~2G{qq) being the dynamically generated quark mass. 



p{k) 



27r2' 



(3) 



Now we consider a quark matter confined in a given sphere that is embedded in the hadronic vacuum and do not 
impose the pressure balance condition between interior of the sphere and external vacuum. We would rather fix 
the radius of the sphere and examine the phase diagram. To incorporate finite-size effects into the thermodynamic 
potential, we replace the density of states p{k) with the MRE density of states PMKE{kj m, R), 



PMRE{k,m,R) 



27r2 



1 



6n^ 



fs 



JkR)' 



where R is the radius of the sphere and the functions. 
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(5) 
(6) 



correspond respectively to the surface and the curvature contributions to the fermionic density of states. The param- 
eter m in Eq. Q should be identified with the mass of quarks in the interior of the sphere. It should be also noted 
that the precise expression of the function fc{x) has not been derived within the MRE framework. Here we used the 
assumption proposed in Ref . ^] . We show the momentum dependence of the MRE density of states in Fig. ^ One 
can see that the MRE density of states is suppressed by the finite-size effects and the suppression is more pronounced 
at low momenta. It can also be seen that the MRE density of states takes unphysical negative values at low momenta 
because of the neglect of the higher order terms in 1/R. Thus, we shall introduce the infrared cutoff Air that is 
defined by the largest zero of pmke,{^ir^ m, R), as is shown in Fig. 

Thus, within the MRE approximation, the thermodynamic potential including the finite-size effects is expressed as 
follows. 



4G 



dkp 



MRE 



IR 



E, 



f3 



1 



(7) 



To obtain the phase diagram we minimize Eq. Q with respect to m. When the value of m at the minimum is 
zero, the system is in the chiral symmetric phase (QGP) and if the value of m is not zero, the system is in the broken 
symmetry phase (hadronic). We then find the phase diagram given in Fig. 13 The outermost curve corresponds to 
R = oo (thermodynamical limit). In this limit, of course our approach recovers the known results[35j. The order of the 
chiral phase transition is second order at small chemical potential and as the chemical potential increases, the order 
of the phase transition is changed form second to first at the tricritical point. In this curve the parts corresponding 
to the second order and the first order transitions are indicated by the dotted and solid lines, respectively. 

To see the effect of finite size in this approach, we show the critical temperature lines for R = 100, 50, and 30 fm 
(from outside to inside). One can see that the critical temperatures are lowered as R is decreased. This behavior is 
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FIG. 1: The momentum dependence of the MRE density of states. 
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FIG. 2: The QCD phase diagram calculated using the MRE density of states given in Eq. (3). The solid and dotted lines 
denote the first order and the second order phase transition, respectively. Each line corresponds to i? = oo, 100, 50 and 30 fm 
from the most outside one, respectively. 



consistent with the enhancement of large fluctuations induced by the finite-size effects. However, we note that the 
tricritical point which exists at i? = oo disappears for all finite values of R. At any value of chemical potential, the 
order of phase transition becomes the first one and the second order phase transition disappears. This is because, as 
we will see below, the MRE cannot describe the second order phase transition. 

These resiilts obtained by the simple application of the MRE to the thermodynamic potential arc physically un- 
reasonable for the following reasons. First, the reduction of the critical temperature seems to be too strong. Second, 
more critical one from physical point of view, is the disappearance of the second order phase transition. In general, 
large fluctuations of the order parameter could increase the order of the phase transition, for instance, from first to 



4 



R=infinity 
R=20 fm 
R=10fm 



< 
3 



-0.15 



-0.16 I 1 1 1 1 

0.2 0.4 0.6 0.8 

m/Ayv 

FIG. 3: The R dependence of the thermodynamic potential. The gradient at m = is positive when the finite-size effects are 
taken into account. 



second. Therefore, it is unnatural that the second order phase transition was converted to the first order one. This 
unphysical behavior indicates the necessity of some modification of the present simple application of the MRE to the 
NJL model. 

To clarify the problem mentioned above, let us take a closer look on the self-consistency condition (SCC) including 
the MRE density of states. The derivative of lo with respect to m is given by 



did 

dm 
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(8) 



[the term proportional to dAm/dm vanishes by virtue of pmre(Air, to, R) — 0]. For the second order phase transition 
to exist, it is necessary that this derivative vanishes at m = 0. However, as has been pointed out in Ref. [sH. m — 
is not the solution to the SCC because of 



m^a (Jm 



(9) 



so that 



duj 
dm 



> 



(10) 



m— 



except for R = oo. That is, the thermodynamic potential has always positive gradient at to = 0, hence the thermo- 
dynamic potential (|7|) cannot describe the second order phase transition. In Fig. |31 the thermodynamic potential at 
(T, /i) = (100, 0) MeV is plotted as a function of to for i? = oo, 20 and 10 fm. 



III. MODIFICATION OF THE DENSITY OF STATES 



In the proceeding section, we examined the T-fi phase diagram, looking for the minimum of the thermodynamic 
potential. This sort of the application of the MRE density of states has been employed in Refs. I^H^- The second 
order phase transition suddenly disappears soon after taking the finite-size effect into account. Thus, the phase 
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diagram with finite size does not converge to that with infinite size even if the radius R is increased. It follows that 
we cannot take the thermodynamic limit in the MRE approximation. Moreover, fluctuations due to the finite-size 
effects is expected to change the order of phase transition from first to second. Hence, it is unacceptable that the 
finite size effect is taken account of the MRE approximation employed in the proceeding section. In this section, we 
reconsider the application of the MRE to resolve this discrepancy. 

The reason why the second order chiral phase transition disappeared in the previous section is the nonexistence of 
the trivial solution to the SCC. The disappearance is due to the mass dependence contained in the MRE density of 
states which prevents the disappearance of the second order phase transition is due to the mass dependence contained 
in the MRE density of states. As is shown in the last of the proceeding section, the SCC does not have the trivial 
solution because of the mass dependence. 

However, there can be another scenario for the application of the MRE. According to the paper by Bailin and 
Blochj23|, the parameter contained in the MRE density of states is not the mass of particle inside the cavities but a 
parameter which characterizes the penetration of the wave function of the quark matter confined in a finite volume 
into the hadronic environment. In the following, we impose specific boundary conditions on quarks and study the 
cases in which the MRE density of states does not depend on the dynamical quark mass. In this scenario, the mass 
parameter contained in the MRE density of states is replaced by the parameter k which is generally independent of 
the dynamical quark mass. Then, the SCC is given by 



7TL f UV jj^ 

— = vl dfcpMRE(fc,K,i?) — 



1 

+ 



(11) 



We do not have the derivative term of the MRE density of states and the SCC obviously has the trivial solution. 
For instance, let us consider the Dirichlet boundary condition by setting k = oo|27j|. Then, we obtain 

hm fs{k/K) - (12) 

ft— ►oo StT 

hm fc{k/K) - -i^, (13) 

and the infrared cutoff is given by 



= m 

In Fig. 01 we show the corresponding phase diagram. The solid and dotted lines represents the critical lines of 
the first order and second order phase transitions, respectively. These lines represent the critical lines of (a) R = oo, 
(b) 20, (c) 10 and (d) 5.7 fm from the out-most one, respectively. Compared to the proceeding result, the finite-size 
effects are suppressed. 

The dashed line represents the trajectory of the tricritical point as a function of R. One can easily see that the 
temperature and the chemical potential of the tricritical point are reduced as the R is decreased. Finally, the tricritical 
point disappears at i? = 12 fm and hence the region of the first order phase transition vanishes. It should be noted 
that the phase diagram calculated with the MRE density of states converges to that in the infinite size in the limit 
of R — > oo, and is consistent with the fact of the increase of fluctuations due to the finite size effects. Thus, this 
approach seems to be more plausible than the proceeding approximation. 

Note that even the second order phase transition disappears at the radius smaller than R = 5.7 fm. This is because 
of the existence of the infrared cutoff and wc will find the phase transition at smaller R by taking account of the 
higher order terms of the MRE density of state. 

If this boundary condition is applicable for the real physical situation created in relativistic heavy ion collisions, 
then the equation of state of the matter changes with the expansion of the system. At the early stage of creation of 
the hot matter at the center, the system size is the order of 10 fm and the critical temperature at = is around 
130 MeV, but after the system expands and reaches the hadronic freezeout stage, then the system size is much larger 
and the critical temperature almost recovers the thermodynamical limit. Therefore, the hydrodynamical description 
of such a situation requires an equation of state which reflects the above finite-size effect. 

The physical environment consistent with the MIT bag model is considered to be realized by setting the Neumann 
boundary condition, k = 0, 

lim/s(fc/«) = 0, (15) 
lirnfcik/n) = (16) 
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FIG. 4: The QCD phase diagram calculated using the MRE density of states given in Eq. (11) with n = oo. The solid and 
dotted lines represents the first order and second order phase transition, respectively. Each line corresponds to (a) R — oo, (b) 
20, (c) 10 and (d) 6 fm from the most outside one, respectively. The dashed represents the trajectory of the tricritical point as 
a function of R. At the radius smaller than R = 5.7 fm, our calculation does not show the chiral symmetry breaking. 



and the infrared cut off is 

In Fig. [SI we show the corresponding phase diagram. The lines correspond to the critical lines for the cases of 
(a) R = oo, (b) 20, (c) 10, (d) 2 and (e) 0.9 fm from the out-most one, respectively. The finite-size effects are more 
suppressed, because the surface term vanishes in this boundary condition, and hence the difference of the critical 
lines between R = oo and i? = 20 fm in Fig. 0]is negligible. The remarkable reduction of the critical temperature is 
observed at smaller size than R = 2 fm. The trajectory of the tricritical point denoted by the dashed line changes 
steeper than the case of the Dirichlet boundary condition. In this boundary condition, the tricritical point vanishes 
at i? = 1.2 fm and the second order phase transition disappears at the radius smaller than R = 0.9 fm. 

Contrary to the previous case, results from this boundary condition indicate that the finite size system with i? = 20 
fm is large enough to take the thermodynamic limit. In this case, it will be reasonable to apply the usual equation 
of state which is given by the thermodynamical limit in hydrodynamical models to describe the heavy-ion collision 
processes. 

So far, we have discussed the behavior of the thermodynamic potential. However, it is not easy to see how the 
second order phase transition is affected by the finite-size effects. For the purpose of looking at characteristic changes 
of the second order transition, we calculate the specific heat defined by 



(18) 
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FIG. 5: The QCD phase diagram calculated using the MRE density of states given in Eq. (11) with n — 0. The solid and 
dotted lines represents the first order and second order phase transition, respectively. Each line corresponds to (a) oo, (b) 20, 
(c) 10, (d) 2 and (e) 0.9 fm from the most outside one, respectively. The three lines corresponding to J? = oo, 20 and 10 fm 
almost degenerate. The dashed represents the trajectory of the tricritical point as a function of R. At the radius smaller than 
R = 0.9 fm, our calculation does not show the chiral symmetry breaking. 



Note that the specific heat shows discontinuity at the critical point of the second order phase transition. 

Figures 6 and 7 show the specific heat of the systems for the cases of k = and k — oo, respectively. Here, we 
chose n = 0.1 GeV for a common quark chemical potential and examined the cases oi R — 100 fm and i? = 10 fm. 
In both cases, the specific heats jump at the corresponding critical temperatures. The jump of the specific heat still 
exists after including the finite-size effects. The result indicates that the second order phase transitions remain second 
order. However, the gap of the specific heat is decreased by the finite-size effects. 

IV. SUMMARY AND CONCLUDING REMARKS 

In this paper, wee have discussed the finite-size effects on the chiral phase transition in QCD in the context of 
the two-fiavor Nambu-Jona-Lasinio model in the mean-field approximation. To incorporate finite-size effects into the 
thermodynamic potential, we employed the multiple refiection expansion (MRE), where finite-size effects are included 
by the modification of the density of state. 

We applied the MRE approximation in different two ways. One is that the "mass" parameter contained in the 
MRE density of state is identified as the dynamical quark mass as is done in l^^l • In this case it is shown that 
only the first order phase transition can take place for any finite size and the critical point disappears. 

The other is that we identify the "mass" parameter as the inverse of the logarithmic derivative of the wave function 
at the surface, following the discussion by Balian and Bloch'27|. This corresponds to the situation that our system 
is confined in a finite volume by an external potential and the boundary condition is fixed independently from the 
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FIG. 6: The specific heat C as a function of T for the cases of i? = 100 fm (•) and i? = 10 fm (□). The mass parameter is 

taken to be « = 0. 
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FIG. 7: The specific heat C as a function of T for the cases of R = 100 fm (•) and J? = 10 fm (□). The mass parameter is 
taken to be /t = oo. 



order parameter of the system.. Then, the order of chiral phase transition changes from the second to the first as the 
chemical potential increases. The temperature and the chemical potential of at the c:ritical point become smaller as the 
system size decreases. This behavior is consistent with the fact that fluctuations is enhanced by the finite-size effects. 
However, quantitatively, we found that the two extreme boundary conditions lead to substantially different results. 
For the Dirichlet boundary condition, the finite-size effect is still large. If this is the case, the critical temperature 
changes from Tc ~ 130 MeV to Tc ~ 180 MeV at zero chemical potential if the size of the system expands from 10 fm 
to 50 fm. Such effects should reflect on the hydro dynamical collective observables, i.e., the V2- On the other hand. 
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for the Neumann boundary condition, the finite-size effects are not so visible and the equation of state is almost the 
same as that obtained from the thermodynamical limit. 

It should be noted that the MRE calculation will lose its validity when the size of our system is very small since 
we have ignored the higher order terms of the MRE approximation. Moreover, the finite-size effects considered in the 
MRE approximation is only the reduction of the density of state. However, in general, we can expect other influences 
due to the finite-size effects; for example, the energy eigenvalue of particles will be also affected. 

We calculated the thermodynamic potential in the mean-field approximation, although it is known that large 
fluctuations near phase transitions invalidate the mean-field approximation. However, the QCD critical temperature 
in the mean-field approximation is qualitatively coincide with that of the lattice QCD calculations. Thus, the results 
obtained in this work would be still qualitatively reliable, unless the fluctuations are extraordinary enhanced by the 
flnite-sizc effects. 

Finally, we comment on the outlook for future studies. In this paper, we studied the static behavior of chiral 
symmetry. However, the finite-size effects would also affect the dynamical behaviors. For instance, the life time of 
fluctuations near phase transitions increases because of the critical slowing down (CSD)[2El2^- The CSD is due to 
the large fluctuations, and hence the finite-size effects will enhance the CSD. 
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